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ABSTRACT 

By exploiting the data of three large surveys (WINGS, NFPS and 
SDSS), we present here a comparative analysis of the Fundamental Plane 
Log(-Re) = Log(cr) + 6 {fi)e + c (FP hereafter) of the early-type galaxies 
(ETGs) belonging to 59 galaxy clusters in the redshift range 0.04 < z < 0.07. 

We show that the variances of the distributions of the FP coefficients derived 
for the clusters in our sample are just marginally consistent with the hypothesis 
of universality of the FP relation. By investigating the origin of such remarkable 
variances we find that, besides a couple of obvious factors, such as the adopted 
fitting technique and the method used to measure the photometric and kine- 
matic variables, the coefficients of the FP are also infiuenced by the distribution 
of photometric/kinematic properties of galaxies in the particular sample under 
analysis. In particular, some indication is found that the FP coefficients intrin- 
sically depend on the particular luminosity range of the sample, suggesting that 
bright and faint ETGs could have systematically different FPs. We speculate 
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that the FP is actually a bent surface, which is approximated by different planes 
when different selection criteria, either chosen or induced by observations, are 
acting to define galaxies samples. 

We also find strong correlations between the FP coefficients and the local 
cluster environment (cluster-centric distance and local density), while the corre- 
lations with the galaxy properties are less marked (Sersic index), weak (color) or 
even absent (flattening). Furthermore, the FP coefficients appear to be poorly 
correlated with the global properties of clusters, such as richness, virial radius, 
velocity dispersion, optical and X-ray luminosity. 

The relation between luminosity and mass of our galaxies, computed by tack- 
ing into account the deviations from the r^/"^ light profiles (Sersic profiles), in- 
dicates that, for a given mass, the greater the light concentration (high Sersic 
index n) the higher the luminosity, while, for a given luminosity, the lower the 
light concentration, the greater the mass. Moreover, the relation between mass- 
to-light ratio and mass for our galaxy sample (with Sersic profile fitting) turns 
out to be steeper and broader than that obtained for the Coma cluster sample 
with r^/^ profile fitting. This broadness, together with the bending we suspect 
to be present in the FP, might partly reconcile the phenomenology of the scaling 
relations of ETGs with the expectations from the ACDM cosmology. 

The present analysis indicates that the claimed universality of the FP of 
ETGs in clusters is still far from being proven and that systematic biases might 
affect the conclusions found in the literature about the luminosity evolution of 
ETGs, since datasets at different redshifts and with different distributions of the 
photometric/kinematic properties of galaxies are compared with each other. 

Subject headings: galaxies: clusters: general — galaxies: fundamental parameters 
— galaxies: structure — galaxies: ellipticals and lenticulars 



INTRODUCTION 



The survey WINGS (IFasano et al. 1120061 ) is providing a huge amount of spectroscopic 
and photometric (multi-band) data for several thousands galaxies in a complete sample of 
X-Ray selected clusters in the local Universe (0.04<z<0.07). Among the other things, line 
indices and equivalent widths (including Mg2 line-strengths) of galaxies are going to be 
available for ~6,000 galaxies, while, for ~40,000 galaxies, we already have at our disposal 
the structural parameters {Re , (fj,)e. and Sersic index n) derived using the automatic surface 
photometry tool GASPHOT (jPignatelli et al. 1120061 ). This put us in a privileged position 
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to analyse the scaling relations of nearby cluster galaxies with unprecedented statistical 
robustness. In this paper we will focus on the Fundamental Plane of early-type galaxies. 



Since its discovery, the FP relation: Log(i?e) = Log(cr) + b (/i)e + c (jPressler et al 



19871 : iDjorgovski fc Davies 1119871 ) has been widely used as a tool to investigate the properties 
of ETGs, to derive cluster dista nces and galaxy peculia r motions (see e. q. th e EN ACS cluster 
survey of 
project of 



<!atgert et al. 



1996 



the SMAC survey of iHudson et al. 1 120011 . and the EFAR 



Wegner et al. 1119961). to per form cosmological tests and compute cosmological 



parameters (see e.g. iMoles et al. I Il998l ) , and as a diagnostic too l of galaxy evolution an d 



M/L variations with redshift (see e.g. iKjaeaegaard et al. I Il993l and IZiegler et al. Ill999l ). 



Most analyses in the literature are based on the comparison between t he FP of distant 



clusters and that of nearby clusters, usually set on the Coma cluster ( |j0rgensen et al 



19961 ). the only one with extensive, homogeneous photometric and spectroscopic data for a 
large sample of ETGs. 

Even if the universality of the FP relation at low redshift has never been actually proven, 
it has been recently claimed that the FP coefficient c0 changes syste matically at increasing 



redshift, from ~ 1.2 at redshift zero to ~ 0.8 at 2; ~ 0. 8 1.3 (Idi Serego et al. 1 12005 



J0rgensen et al. 2006). This change, already predicted by Pahre et al. ( 1998al ). has been 
attributed to the evolution of ETGs with redshift. 

However, the situation is far from being clear, since t he data required to a ssess the 
universality of the FP are still lacking. The SDSS survey (IBernardi et al. II2003I ) first at- 
tempted to face this problem adopting the correct strategy, which must necessarily rest on 
the availability of large galaxy samples. The results of this analysis indicate that the FP is 
a robust relation valid for all ETGs (above the magnitude limit of the SDSS), but its co- 
efficients could depend on the number density of the galaxy environment: the luminosities, 
sizes, and velocity dispersions of the ETGs seem to increase slightly as the local density 
increases, while the average surface brig htnesses decrease. Howey e r, evidences supporti ng 



diffe rent conclusions h ave b een found by Ide la Rosa et al. 
and iKochanek et al. I J2000h . 



fl200l[ ). IPahre et al. I fll998al Jbl) 



In addition, it is still unclear whether ETGs in clusters at the same redshift share the 
same FP, or instead the FP coefficients systematically change as a function of the global 
properties of the host cluster (richness, optical and X-ray luminosity, velocity dispersions, 
concentration, subclustering, etc.). 



^This coefficient is related to the tilt of the FP, represented by the difference 2 — a, that is the deviation 
from the Virial expectation value a — 2. 
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Today, thanks to t he huge observational effort done by wide fi eld surveys, such as SPSS 



flBernardi et al. Il2003h . NFPS fISmith et al. Il2004f ) and WINGS flFasano et al. Il2006f ). the 



study of the FP can be extended to a much larger sample of nearby clusters. Besides the 
data from the SDSS survey, we can now use those from two more surveys (WINGS and 
NFPS) suitably designed to study the properties of nearby clusters. Here we exploit these 
datasets to check whether, at least in the local Universe, the hypothesis of universality of 
the FP turns out to be supported by the observations or not. 

The paper is structured as follows: in Sec. [2] we present our data sample, discussing its 
properties, its statistical completeness and the intrinsic uncertainties associated to the mea- 
sured structural (effective radius), photometric (effective surface brightness) and dynamical 
(central velocity dispersion) quantities involved in the FP relation. In Sec. [3] we present the 
FP for the whole dataset and those of each individual cluster. In Sec. HJ also by means 
of extensive simulations, we investigate the origin of the large spread observed in the FP 
coefficients, showing that the scatter is hardly attributable just to the statistical uncertainty 
arising from the limited number of ETGs in each cluster. In Sec. [5] we explore the behaviour 
of the FP coefficients at varying some galaxy properties (Sersic index, color, flattening), the 
local environment (cluster-centric distance and local galaxy density) and the global proper- 
ties of the host clusters (density, central velocity dispersion, optical and X-ray luminosity). 
Finally, in Sec. [6], we discuss the relations involving the mass and the mass-to-light ratio 
of ETGs in nearby clusters, which are closely linked to the FP, also providing a tool to 
investigate the galaxy formation and evolution. Conclusions are drawn in Sec. [71 Hereafter 
in this paper we adopt the standard cosmological parameters Hq = 70, Qx = 0.7, Qb = 0.3. 



2. THE GALAXY SAMPLE 

The initial galaxy sample has been extracted from 59 clusters belonging to the survey 
WINGS (W). It includes galaxies having velocity dispersion measurements and 'early- type' 
classifications from the surveys SDSS (S) and/or NFPS (N). Effective radius a nd surface 



brightness of galaxies have been measured by GASPHOT (jPignatelli et al. II2006I ). the soft- 
ware purposely devised to perform the surface photometry of galaxies with threshold isopho- 
tal area (at 2xrmsbkg) larger than 200 pixels in the WINGS survey (Pignatelli et al. in 
preparation). The central velocity dispersions have been extracted from the catalogs pub- 
lished by the surveys NFPS (52 clusters in common with WINGS) and SDSS (14 clusters 
in common with WINGS). The clusters in common between NFPS, SDSS and WINGS are: 
A0085, A119, A160, A602, A957x, A2124, and A2399. 



A careful check of morphologies, performed both visually and using the automatic tool 
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MORPHOT (Fasano et al. in preparation; again purposely devised for the WINGS survey), 
allowed us to identify in both datasets several early-type spirals, erroneously classified as E 
or SO galaxies (~ 8% of the whole sample). Besides these, we also decided to exclude from 
the present analysis the galaxies with central velocity dispersion a < 95 km s~^(see Sec J2.2l) 
or total luminosity My > —18. The final sample sizes are: iVvi/+Ar=1368; Nw+s='^82; 
Nw+N+s=^550 (100 objects in common between W+S and W+N). The median number of 
ETGs per cluster is Nmed = 23. For each cluster. Table [1] reports the number of galaxies in 
the two samples (W+N and W+S; columns 8 and 9, respectively) and that of galaxies in 
common (W+[N&S]; column 10). 

The table also reports some salient cluster properties: average redshift (column 2; from 
NED), velocity dispersion (S) of galaxies around the average red shift (column 3; again frona 



NED ), X-ray (0.1-2.4 keV) luminosity in ergs s ^ (column 4; from lEbeling et al. Ill996l . Il998 



20001 ) ■ total absolute magnitude in the V-band (column 5; from the W INGS deep catalogs). 



radius -R200 in Mpc (column 6; from E, following iPoggianti et al. II2OO6I ) and absolute V-band 



magnitude of the brightest cluster member (column 7; again from the WINGS catalogs). 

It is worth stressing that, even though our sample of ETGs is the most sizeable among 
those used till now to study the FP of nearby clusters, it is still far from being complete from 
a statistical point of view. In particular: (i) the surface photometry is available just for the 
galaxies in the region of ~ 35 x 35 arcminutes around the cluster center (the regions mapped 
by the CCD images of the WINGS survey); (ii) the SDSS and NFPS surveys have provided 
velocity dispersions just for subsamples of the WINGS ETGs, each survey according to the 
proper selection criteria (see Sec j2.2p : (iii) a couple of clusters with SDSS velocity dispersions 
are just partially mapped by the survey. 



2.1. The WINGS photometry 

The WINGS survey has produced catalogs of deep photometry and surface photometry 
for 77 nearby clusters. For several thousands galaxies per cluster the deep catalogs contain 
many geometrical and aperture p hotometry data (Varela et a l. 2008, A&A, in press.), derived 



by means of SExtractor analysis (IBertin and Arnouts Ill996l ). The surface photometry cata- 



logs contain data for several hundreds galaxies per cluster (those with isophotal area greater 
than 200 pixels) and have been produced by using the previously mentioned tool GASPHOT. 
For each galaxy it performs seeing convolved, simultaneous, Sersic law fitting of the major 
and minor axis growth profiles, thus providing Sersic index n, effective radius Re and average 
surface brightness (/i)e , total luminosity, flattening and local sky background. The data and 
the associated uncertainties are discussed in Pignatelli et al. (2008, in preparation). The 
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average quoted rms uncertainties of Re and (/i)e are ~ 15% and ~ 10%, respectively. The 



surface brightnesses have been corrected for galactic extinction (jSchlegel et al. Ill998l ) and 
cosmological dimming (using the average redshifts of the clusters), while the K-corrections 
have not been considered. Effective radii have been transformed from arcseconds to Kpcs 
using the cosmological parameters given in Section [TJ 

It is worth stressing that just a few dozens of galaxies per clusters, out of the several 
hundreds for which WINGS provides surface photometry parameters, can be included in the 
final sample, due to the morphological constraint (early-type) and the cross matching with 
the available velocity dispersion data. 

In Figure [1] we compare total magnitudes and effective radii derived by GASPHOT 
(Sersic's law fitting) with the corresponding quantities derived by the SDSS surface pho- 
tometry using the de Vaucouleur's r^^^ law. The 407 ETGs in common between the SDSS 
and WINGS surveys (including galaxies with a < 95 or My > —18) are shared among 14 
different clusters. The figure clearly illustrates how the surface photometry parameters are 
strongly influenced by the adopted fitting procedure. In particular, in our case, the strong 
dependence of both ALog(Re) and AV on the Sersic index n is largely expected due to the 
different amount of light gathered in the outer luminosity profiles by the R^^^ and Sersic law 
extrapolations. However, it is worth noting in Figure [1] that, even for ?t,=4 (dotted lines in 
the figure) the GASPHOT and SDSS surface photometries give different results, the last one 
producing slightly fainter and smaller galaxies. To this concern, according to the SDSS-DR6 
documentation, both the effective radii and the total luminosities provided by SDSS for 
galaxies in crowded fields (as the clusters are) turn out to be more and more underestimated 
at increasing the galaxy luminosity. In the magnitude range typical of our galaxy sample 
(~15<V<~18) we expect these biases to be of the order of -0.05 and 0.05 for ALog(Re) 
and AV, respectively. While for ALog(Re) the expected bias could be enough in order to 
explain the discrepancy in the figure (upper panel), for AV (lower panel) it would be largely 
insufficient. The residual discrepancy (AV~0.15) is likely attributable to the difference be- 
tween the fitting algorithms used by SDSS (2D - pixel by pixel) and GASPHOT (major and 
minor axis growth profiles; see Pignatelli et al. 2006 for a discussion of the advantages of 
this fitting procedure). 



2.2. The kinematical data 

The central velocity dispersions a of the ETGs have been taken from the published data 
of the NFPS and SDSS-DR6 surveys. It follows that the completeness is strongly affected 
by the selection criteria adopted in these surveys. In particular, the SDSS survey defines 
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Fig. 1. — (Upper panel): difference between the effective radii derived by the SDSS surface 
photometry using the de Vaucouleur's r^/^ law and those derived by GASPHOT using the 
Sersic's law, as a function of the Sersic index n (by GASPHOT), for the 407 ETGs for which 
both the SDSS and the WINGS surveys provide surface photometry parameters. (Lower 
panel): as in the upper panel, but for the total V-band magnitudes. In this case, the 
SDSS V-ban d magnitudes are obtain ed from the r'-band ones using the conversion formula 



proposed by lFukugita &: et al. Ill996 
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ETGs those objects having both a concentration index Rgo/R^o > 2.5 (in the i* band) and 
a very good r^/^ de Vaucouleurs hght profile, while the ETGs of the NFPS survey have 
been selected on the basis of their colors, using a narrow strip around the color-magnitude 
diagram. Both criteria might lead to exclude from the samples the brightest cluster galaxies 
(BCGs), which are actually laking in the SDSS sample. Moreover, in the SDSS survey the 
velocity dispersions are measured only for spectra with signal-to-noise ratio S/N > 10 (high 
average surface brightness) and some clusters are not fully mapped by the survey strips. 
We will see that such different selection criteria produce systematic differences in the FP 
coefficients derived for the two samples. 

It is worth pointing out that in the originally submitted version of this paper (arXiv0804 . 1892D) 
we used SDSS velocity dispersion data from a previous release of the survey (SDSS-DR4) 
and that the differences between the velocity dispersions given in DR4 and DR6 are not 
negligible, especially for small values of a (see Figure [2]). This is the reason why many fig- 
ures and tables, as well as some findings we report here (mainly concerning the difference 
between the FP coefficients of the SDSS and NFPS samples) are slightly different from the 
corresponding ones reported in the previous version of the paper. Still, we decided to keep 
that version unchanged on the babbage (just slightly modifying the title) in order to show 
how much a correct determination of the physical quantities involved in the FP (especially 
a) is critical in drawing any conclusion from the FP tool. 

All the available velocit y dispersions have been homogenized to the uniform aperture 
Re/8, following the recipe of I J0rgensen et al~ ( 1995 ). The estimated uncertainty for both 



surveys is in the range 7 -r 10%. 

In Figure [3] we plot the difference Log((jAr)-Log(crs') versus Log((7Ar) for the 100 galaxies 
of our sample in common between the NFPS and SDSS samples. The rms scatter of the 
Log(crAr) vs Log(cr5) relation is ~ 0.05, equivalent to an uncertainty of ~ 12% in the common 
velocity dispersions. Again there is a systematic deviation between the two datasets at low 
velocity dispersions (cr < 95 km s"^). 

In the following, to avoid any possible bias in the comparison of the FP of clusters, we 
have excluded from our analysis the objects with a < 95 km s~^. Moreover, when dealing 
with the global (W-f-N-l-S) galaxy sample, the average velocity dispersion a = {aN + crs)/2 
have been assigned to the galaxies in common between NFPS and SDSS. 
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Log((7DR6) 



Fig. 2. — Comparison between Log{aDR6) and Log^aDRi) for 523 SDSS galaxies originally 
selected in the fields of our WINGS survey. Note the systematic offset between the DR4 and 
DR6 releases at low velocity dispersions. 
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Fig. 3. — The difference Log((T5)-Log((j]v) versus Log((T5) for the 100 ETGs of our sample 
in common between the NFPS and SDSS surveys. Note the systematic offset at low velocity 
dispersions, which led us to restrict our sample to galaxies with cr > 95 km . 
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3. FITTING THE FP 

It is well known that the values of the FP coefScients vary systematica lly at varying 



the adopted fitting algorithm fIStrauss &: Willick Ill995l : iBlakeslee et al. II2002I ) and that the 



choice of the algorithm actually depends on the particular issue under investigation (relation 
among the physical quantities, linear regression for distance determination, etc.). Here we 
tried two di fferent algorithms to get the best fit of the FP: 1) the program MIST, kindly 



provided by iLa Barbera et al. I (120001 ). which is a bisector least square fit, coupled with a 
bootstrap analysis providing a statistical estimate of the errors of the FP coefficients; 2) a 
standard fit minimizing the weighted sum of the orthogonal distances (ORTH hereafter). 
Both algorithms account, in different ways, for the measurement errors on the variables 
Log(i?e), (Ai)e, and Log((T). MIST considers an average covariance matrix that includes the 
variances of the errors in all parameters and their mutual correlations (such as rmsLog(Re) vs. 
rms ). On the other hand, ORTH takes into account the errors of individual measures 
in a standard analysis. In Table [2] we report the FP coefficients derived from the two 
fitting algorithms for the global galaxy sample (first two lines) and for the NFPS and SDSS 
samples separately (lines 3-4 and 5-6, respectively). In the same table (lines 8-9) we report 
for comparison the FP coefficients obtained with both MIST and ORTH fitting algorithms 
for a sample of 80 E TGs in the Coma cluster (photometric and kinematical data from 
J0rgensen et al. lll995l [JORG]l The column labeled with A''^ in the table reports the number 



of galaxies used in each fit. 

Besides the best fitting algorithm, the FP coefficients might also be systematically 
infiuenced by the technique adopted to measure the effective radius and surface brightness 
of galaxies (1D/2D light profile fitting with de Vaucouleurs/Sersic laws). Lines 5 and 7 of 
Table [2] report the MIST FP coefficients obtained for the galaxy sample in common between 
WINGS and SDSS, using alternatively the two surface photometry data sets (see in Figured] 
the comparison among them and in Section 12.11 the description of the WINGS and SDSS 
surface photometry techniques). It is evident that, at least in our case, the infiuence of the 
adopted surface photometry technique on the FP coefficients turns out to be negligible. 

Table [2] shows that different fitting algorithms (and, possibly, surface photometry tech- 
niques) lead to somewhat systematic differences in the FP coefficients. In particular, the 
values of a obtained using the MIST fit are in general slightly smaller than those coming 
from the orthogonal fit. This means that, in order to perform a correct comparison of the FP 
results, it is advisable to adopt homogeneous FP fitting and (perhaps) surface photometry 
techniques. 



However, in the present analysis, we do not focus on the 'true' values of the FP coeffi- 
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cients. Instead, we will concentrate on their possible variation as a function of both galaxy 
and cluster properties. In other words, rather than in obtaining the best possible fit for a 
given application of the FP, we are interested in investigating the FP systematics, once both 
the fitting algorithm and the surface photometry technique have been chosen. Hereafter we 
adopt the MIST bisector fitting algorithm and the WINGS-GASPHOT surface photometry. 
The last choice will allow us to account for the structural non- homology of galaxies (Sersic 
index, see Section [6]), while the former one will provide FP coefficients useful for distance 
determination of farther clusters. However, using the ORTH fitting algorithm, we will also 
provide in Section 15.31 a recipe for the V-band FP, useful to define the physical relation 
among the quantities involved in it. 

Comparing each one another the FP coefficients given in Table [21 we easily realize that, 
besides the obvious dependence on the fitting algorithms, a further dependence exists on the 
galaxy sample, even adopting the same fitting algorithm (MIST) and surface photometry 
technique (WINGS-GASPHOT). In particular, the a coefficient, which is related to the so 
called 'tilt' of the FP, is noticeably different for the three data samples, even if the rms 
scatter in Log(i?e) is always ~ 0.05 (which implie s an u ncertainty of ~ 12%), a value just a 



bit larger than that reported in I J0rgensen et al. I (119951 ) (~ 11%). 



In panel (a) of Figure H] we show the MIST bisector fit of the FP obtained for the whole 
W+N+S dataset (see line 1 of Table [2]) using two different colors for the W+N and W+S 
data samples (respectively black and grey; green in the electronic form). Note the cut shown 
by SDSS data at large values of hogi^Re), which is obviously due to the bright end cut of the 
survey. 

In Figure [5] we plot the FP of the individual clusters, again using for reference the 
coefficients derived from the fit of the whole W+N+S data sample. Columns 3-8 of Table [3] 
report the best fit coefficients (and the associated uncertainties) of each cluster, obtained 
with the MIST algorithm. Even from a quick look of both Figure [5] and Table [3l it is clear 
that the global fit does not seem to be a valid solution for all clusters. 

The average values (with their uncertainties), the standard deviations and the median 
values of the (MIST) FP coefficients of the clusters in the global sample and in the samples 
W+N and W+S are reported in Table. |H From this table we note that: (i) even if the 
scatter is large, the average values of the FP coefficients appear systematically different 
(well beyond the expected uncertainties) in the W+N and W+S cluster samples, confirming 
the dichotomy already noted in Table [H (ii) when just clusters with A^^ > A^med(= 23) are 
considered, the standard deviations of the distributions of the FP coefficients decrease only 
slightly, suggesting that the large scatter cannot be ascribed to the statistical uncertainties 
related to the (sometimes) small number of ETGs in our clusters. 
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Fig. 4. — (panel a): The FP of the W+N (black dots) and W+S (gray dots; green in the 
electronic version) data samples; (panel b) The W+N+S FP for E (black dots) and SO (gray 
dots; green in the electronic version) galaxies. In both panels we used for reference the FP 
coefficients derived from the best fit of the global dataset (W+N+S). The two-sided arrow 
in panel (a) roughly defines, through the Faber- Jackson (L— a) relation, the direction of 
constant luminosity (or a; see Section l42l) . 
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Fig. 5. — The FPs of individual clusters in our sample are plotted using the MIST best fit 
solution found for the global W+N+S galaxy sample. 
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Fig. 5. — The FPs of individual clusters (continued). 
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4. ORIGIN OF THE SCATTER OF THE FP COEFFICIENTS 

We test two different hypotlieses to explain tlie differences between the W+N and W+S 
samples and, in general, the large observed scatters of the FP coefficients: (i) they are 
simply due to the statistical uncertainties of the fits; (ii) they are artificially produced by 
the different criteria used to select ETGs in the NFPS and SDSS surveys. 



4.1. Consistency with statistical uncertainties 

First we test the 'null hypothesis' that the observed scatter is merely consistent with 
the statistical uncertainties of the fits. To this aim, using all the galaxies in our sample, we 
produced two different sets of simulated clusters. In the first set we generate mock clusters 
with number of galaxies {Ng) progressively increasing from Log{Ng) = 1 to 2 (step 0.1) and 
fit each mock cluster with MIST. Figure [6] shows the average values of the FP coefficients 
and the corresponding standard deviations as a function of Log(A^g). Note that, since for 
each value of A''^ the whole sample of 1550 galaxies is used to randomly extract as many 
mock clusters as possible avoiding galaxy repetitions, the number of mock clusters increases 
at decreasing A''^, thus resulting in almost constant error bars of the average FP coefficients 
and of their variances. 

In the second set of simulations we produced 100 toy surveys, each containing 59 clusters 
obtained by sorting randomly the whole galaxy sample and taking sequentially the same 
number of galaxies per cluster as the real survey (thus avoiding galaxy repetitions)!^. Then, 
using the MIST algorithm, we evaluate the FP coefficients of each mock cluster and, for each 
mock survey, we compute the average and median values of the coefficients, together with 
their standard deviations. Finally, we compare the distributions of the average coefficients 
and their variances in the mock surveys with the corresponding values of the real survey. 

Figures [7| and [8] illustrate the conclusions of the two sets of simulations. The first set 
has been used to compute (with the equations given in the right panels of Figure [6]) the 
error bars in Figure [3 The left panels of this figure report the FP coefficients of our 'real' 
clusters versus the number of galaxies in each cluster, while the histograms on the right side 
of each panel show the corresponding distributions (see figure caption for more details). The 
error bars are used to compute the reduced Chi-square values (reported in the figure; in 
our case z/=58) of the differences between the coefficients of the individual clusters and the 



^Note that, in this way, we imphcitly assume that the probabiHty distributions of photometric/kinematic 
properties of galaxies are the same in all clusters and correspond to those of the global galaxy sample. 
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Fig. 6. — Average values (left panels) and standard deviations (right panels) of the FP 
coefficients as a function of the number of galaxies (Ng) for mock clusters randomly extracted 
from the whole galaxy sample (sec text for details). The dashed lines in the left panels 
correspond to the FP coefficients obtained fitting altogether with MIST the 1550 galaxies in 
our sample. The full lines in the right panels correspond to the simple exponential functions 
we used to compute the standard deviations as a function of Log{Ng) (see the equations in 
each panel). 
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corresponding coefficient of the global galaxy sample (dashed lines in the figure). Apart from 
the a coefficient {P^ ~0.965), they correspond to very high values of the rejection probability 
[Py >0.995) that the coefficients of the individual clusters are randomly extracted from the 
same parent population. 

Figure [8] shows that the average values of the FP coefficients for the clusters of the real 
survey are just marginally consistent with the corresponding distributions obtained with the 
simulations of mock surveys (upper panels), while the distributions of variances are more or 
less in agreement with the real ones (lower panels). 

The two sets of simulations indicate that the observed scatter is not accounted for by 
the statistical uncertainties of the fits and that the real clusters cannot be merely assembled 
by random extraction of galaxies from the global population. The left panels of Figure [7] 
also clearly illustrate the systematic differences between the FP coefficients of the NFPS 
and SDSS samples already quoted in Tables [2] and HI To this concern, the two-sample 
Kolmogorov-Smirnov test, applied to the black and open+gray (green) samples in that Fig- 
ure, provides rejection probabilities of 0.998, 0.530 and 0.986 for the left panels coefficients 
a, h and c, respectively. 



4.2. Dependence on galaxy sampling 

As an early test, we wanted to investigate the hypothesis that the observed scatter 
of differences in the FP coefficients are the result of blending E and SO galaxies, with the 
knowledge that the E/SO ratio varies (for example) as a function of local density. We 
verified that the FP computed separately for the elliptical and SO galaxies are practically 
indistinguishable (all FP coefficients differ by < 3%; see also the right panels of Figure H]). 
This allows us to rule out the hypothesis that the observed scatter is induced by different 
E/SO fractions in the different samples. 

On the other hand, we have seen in Tables [2] and H] that the FP coefficients of the clusters 
in the W+N and W+S data samples are systematically different from each other (see also 
the left panels of Figure [7] and the last sentence of the previous sub-section). It is therefore 
natural asking which is the origin of such systematic difference. 

Since for all galaxies in the sample the surface photometry data come from WINGS+GASPHOT, 
we could be tempted to conclude that the differences we found are due to some systematic off- 
set between velocity dispersion measurements from the NFPS and SDSS surveys. However, 
this possibility is definitely ruled out by Figure [3] (Sec. [2^ . which shows that the agreement 
between the two velocity dispersion surveys is fairly good, at least for a >95km s~^. Indeed, 
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Fig. 7. — (Left panels): FP coefficients of our clusters versus the number of galaxies in 
each cluster. The black and gray (green in the electronic version) dots refer to clusters with 
only NFPS and SDSS galaxies, respectively. The open dots represent the seven clusters in 
common between the NFPS and SDSS surveys. The dashed lines correspond to the FP 
coefficients obtained fitting altogether with MIST the 1550 galaxies in our sample. The 
histograms on the right represent the distributions of FP coefficients in our cluster samples. 
Black, gray (green in the electronic version) and open histograms have the same meanings as 
in the left plots and are cumulated inside each bin; (Right panels): as in the left panels, but 
using only galaxies with My < —19.5. Note that in the left panels (global galaxy sample) 
the NFPS and SDSS clusters have quite different distributions of the coefficients, while in 
the right panels (just galaxies with My < —19.5) the distributions of the two samples are 
consistent among each other. 
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Fig. 8. — Histograms of the average values of the FP coefficients (upper panels) and standard 
deviations (lower panels) for the 100 toy surveys (see text for details). The dashed lines in 
the histograms mark the corresponding values obtained from the real survey (see Table. H]). 
These values are also reported in the panels, together with the probabilities that they are 
randomly extracted from the underlying histograms. Note in the upper panels that the 
average values of the FP coefficients for the clusters of the real survey are just marginally 
consistent with the corresponding distributions obtained from the mock surveys. 
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we have also verified that the FP coefficients of the galaxy sample in common between NFPS 
and SDSS, obtained using alternatively the two velocity dispersion data sets do not differ 
significantly. 

Thus, we are left with the last possibility: that the systematic FP differences between 
the NFPS and SDSS clusters are due to the different distributions of photometric/kinematic 
properties of galaxies in the two samples. The da n ger o f selection biases i n this game 



has already bee n emp hasized by iLynden-Bell et al. I (119881 ) , IScodeggio et al. I (119981 ) , and 



Bernardi et al. I (120031 ). who showed that robust fits of the FP can be obtained only for 
galaxy samples complete in luminosity, volume, cluster area coverage and stellar kinematics. 
The panels (a) and (b) of Figure [9] respectively show the projection of the FP on the surface 
photometry plane ((yu)e — Log(i?e) ; Kormendy relation) and the Color- Magnitude diagrams 
[My — [B — V)] for the NFPS and SDSS surveys. Both figures show that the two galaxy 
samples have quite different distributions of the photometric quantities involved in the FP 
parameters. This is even more evident in the panel (c) of the same figure, where the face-on 
view of the FP of the global sample is shown, together with the loci corresponding to some 
constant values of the quantities involved in the FP (dotted lines). Note that, both in the 
Color- Magnitude and in the Kormendy diagrams, the W+S galaxy sample turns out to be 
(on average) fainter than the W-l-N sample, especially in the small size region. This is likely 
a direct consequence of the rules the two surveys adopt to select early-type galaxies (see 
Section [22D. 

The fact that such differences in the galaxy sampling produce the observed differences in 
the FP coefficients is shown in Figure [TOl In the upper panels of the figure the a coefficient 
of the FP seems to be anti-correlated with the average values of luminosity, radius and 
velocity dispersions of galaxies in the clusters. The same, but (obviously) with positive CCs, 
happens for the coefficient c (not reported in the figure). We see from the lower panels in the 
figure that, if we cut the data samples at higher luminosity. My = —19.5, these correlations 
disappear, since in this case the two data samples are more homogenous. 

This is also conffimed by the right panels of Figure [71 where the plots in the left panels 
are repeated using only galaxies with absolute magnitude My < —19.5. Indeed, the two- 
sample Kolmogorov-Smirnov test, applied to the black and open-|- gray (green) samples in 
that figure, provides rejection probabilities of 0.475, 0.308 and 0.318 for the right panels 
coefficients a, h and c, respectively (compare these values with those given in the last sentence 
of Section 0|) 

A final, quantitative estimate of the dependence of the FP coefficients on the luminosity 
distribution of the galaxy sample is provided by Figure [TH where we report the FP coefficients 
obtained for different values of the faint and bright luminosity cut-off applied to the NFPS 
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Fig. 9. — (Panel a): the (yu)c — Log(i?e) relation for the W+N (black dots) and W+S (gray 
dots; green in the electronic version) samples. The big dark-gray (magenta in the electronic 
version) and white dots represent the average surface brightnesses of the two galaxy samples 
in different bins of Log(i?e); (Panel b): the Color-Magnitude diagrams for the galaxies of 
the W+N and W-l-S samples. Symbols are as in panel (a); (Panel c): face-on view of 
the FP obtained for the global galaxy sample. The dotted straight lines mark the loci 
corresponding to some constant values of luminosity, surface brightness, effective radius and 
velocity dispersion. Symbols as in the previous panels. 
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Fig. 10. — (Upper panels) the FP coefficient a vs. the average values of effective radius 
(< Log{Rc) >), luminosity (< My >) and central velocity dispersion (< Log{a) >) of the 
galaxies in each cluster. Symbols are as in Figure [71 (Lower panels) the same plots, but 
using only galaxies with My < —19.5. 
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and SDSS samples. 

The left panels of Figure [TT] show that, for both the NFPS and the SDSS sample, the 
coefficients a and h decrease at increasing the faint luminosity cut. This effect can be at least 
partially explained by the very geometry of the PP. In fact, in the edge-on representation 
of the PP, any luminosity cut in the galaxy sample translates, through the Paber- Jackson 
(L— 0") relation, in a sort of 'zone of avoidance' delimited by a line of constant luminosity (or 
0"), whose direction is roughly indicated by the two-sided arrow in PigureH] (panel a). This 
Malmquist-like bias reduces the PP slopes along the directions of a and (/i)e for both faint- 
and bright-end luminosity cuts. Pigure [12] illustrates this 'geometrical' eSect. It is similar to 
Pigure \TT\ but compares the W+N+S sample (black dots) with a mock sample of 10,000 toy 
galaxies (grey dots; green in the electronic version of the paper) randomly generated around 
the same (W+N+S) PP, according to the 'true' distributions (and mutual correlations) of 
(/i)e -Re and cr. The right panels of Pigure [11] show that this bias actually works for the 
bright-end luminosity cut just in the case of the NPPS sample. Instead, the PP coefficients 
of the SDSS sample display a rather peculiar behaviour. Por the faint-end cuts they show 
trends similar to those of the NPPS sample, but more pronounced. Instead, they do not 
seem to depend at all on the bright-end cuts (right panels), remaining significantly higher 
than in the case of the NPPS samples over the range of cut-off luminosities. This behaviour 
suggests that, besides the luminosity cut-off, other causes may contribute to tell apart the 
two samples. 

The Pigure [13] helps to clarify this point. It shows the edge-on PP as it appears along the 
direction of luminosity. This particular projection highlights a weak feature of the PP that 
otherwise would be completely masked, suggesting the existence of a sort of warping (black 
curve in the figure; red in the electronic version). Although just hinted in the bright part of 
the luminosity function, this feature looks a bit more evident in its faint-end, which in our 
sample is dominated by SDSS galaxies. To this concern, it is worth noticing that this faint 
luminosity warp can hardly be attributed to a possible upward bia s of the SDSS low velocity 



dispersion measurements, since, according to ISmith et al. I (120041 ) (see also Pigure [3|), such 



a bias should in case work in the opposite direction. The different shape of the NPPS and 
SDSS samples in this particular projection of the PP explain the reasons why: (i) for the 
SDSS sample the coefficient a turns out to be always greater than in the case of the NPPS 
sample; (ii) the faint-end luminosity cut influences the coefficient a of the PP more for the 
SDSS than for the NPPS sample; (iii) the bright-end luminosity cut does not influence the 
PP coefficients of the SDSS sample. 

Although these analyses would benefit from a more robust statistic, they lead us to 
suggest that the PP is likely a curved surface. This fact has been recently claimed by 
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Fig. 11. — The FP coefficients a and b obtained with MIST for different values of the faint 
and bright luminosity cut-off applied to the NFPS galaxy sample (black dots) and to the 
SDSS sample (gray dots; green in the electronic version). 
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Fig. 12. — similar to Figure [TH but comparing the W+N+S sample (black dots) with a 
mock sample of 10,000 toy galaxies (grey dots; green in the electronic version of the paper; 
see text for more datails) 
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Fig. 13. — Edge-on FP as it appears along the direction of luminosity. Symbols are as in 
panel (a) of Figure HI The black curve (red in the electronic version) just represents a naive 
fitting we made in order to enhance the warp-like feature of this particular projection of the 
FP. 
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Desroches et al. I (120011 ) and, in the low-luminosity region, may actua lly indicate a first 



hint of the connection between the FP of gia nt and dwarf ellipticals (iNieto et al. I Il990 



Held et al. I l200ll . iPeterson and Caldwell I Il993l ) . In Section E] we will also present a further 
hint of the existence of the high luminosity warp of the FP suggested by Figure [T31 

It is important to stress that the possible curvature of the FP may give rise to different 
values of its coefficients when different selection criteria, either chosen or induced by observa- 
tions, are acting to define galaxy samples. This fact represents a potentially serious problem 
when the goal is to compare the tilt of the FP at low- and high-redshifts, since it implies 
that a reliable comparison can be done only if galaxy samples at quite different distances 
share the same distributions of the photometric/kinematic properties, which is indeed not 
usually the case. 

Finally, we note that, according to the Chi-square values reported in the right panels of 
Figure [71 the scatter of the FP coefficients is poorly consistent with the expected statistical 
uncertainties, even after having reduced the annoying dichotomy between the NFPS and 
SDSS data samples. This fact suggests that at least part of the observed scatter must be 
somehow 'intrinsic' and resulting from a 'true' dependence of the FP coefficients on the 
galaxy properties and/or on the local environment and/or on the global cluster properties. 
The huge amount of data available from the WINGS photometric catalogs allows us to 
perform for the first time this kind of analysis. 



5. SYSTEMATICS OF THE FP COEFFICIENTS 

In order to reduce the luminosity- driven bias of the FP coefficients illustrated in the 
previous Section 14.21 we decided to use in this section only galaxies with My < —19.5 
(A^g=1477). Even if this luminosity cut-off does not remove completely the systematic FP 
differences arising from the different sampling rules of the NFPS and SDSS surveys (see 
Figure [TT]) . we guess it is able at least to reduce them down to an acceptable level. 



5.1. FP versus galaxy properties and local environment 

In Section IT2] we have already shown that the FP coefficients do depend on the average 
luminosity of the galaxies in the sample and, therefore, on the average values of size and 
velocity dispersion (see Figure [TUl) . These dependences concern the very shape of the FP 
relation, since they involve the physical quantities defining the relation itself. Now, besides 
these 'first order' dependences, we want to check whether the FP relation varies with other 
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galaxy properties or the local environment. In particular, as far as the galaxy properties are 
concerned, we test the {B — V) color, the Sersic index Log(n) and the axial ratio b/a, while 
the cluster-centric distance Dec (normalized to i?2oo and the local density Log(p)0 are 
used as test quantities of the local environment. 

A simple way to perform such kind of analysis is to correlate the tes t quantities with 



the r esiduals of the FP relation obtained for the global galaxy sample (jj0rgensen et al. 



19961 ). For instance, in Figure O the FP residuals are reported as a function of both Dec 
and Log(p). From this figure one would be led to conclude that these two parameters do 
not influence at all the FP coefficients. However, this method would be intrinsically unable 
to detect any correlation if the barycentre of galaxies in the FP parameter space does not 
change at varying the test quantity. In fact, in this case, any change of the slope alone would 
produce a symmetric distribution of the positive and negative residuals, keeping zero their 
average value. For this reason, we preferred to perform the analysis by evaluating the FP of 
galaxies in different bins of the test quantities. Moreover, in order to get similar uncertainties 
of the FP coefficients in the different bins, we decided to set free the bin sizes, fixing the 
number of galaxies in each bin {Nbin). 

In Figure [15] the average values of the the FP coefficients in different bins of the test 
quantities are plotted as a function of the median values of the quantities themselves inside 
the bins. The panels also report the correlation coefficients (CC) of the different pairs of 
bin-averaged quantities. In these plots we set Nhin=i50 and assumed the centers of the 
clusters to coincide with the position of the BCGs. However, the trends and the correlation 
coefficients in the figure remain almost unchanged if we set (for instance) iVbj„=200 and 
assume that the cluster centers coincide with the maximum of the X-Ray emission. 

It is worth stressing that the FPs we obtain with the outlined procedure for each bin of 
the test quantities do not refer to real clusters. They are actually relative to ideal samples 
for which some galaxy /environment property is almost constant (for instance: constant local 
density) . 

At variance with the conclusions one could draw from Figure [HI it is clearly show in 
Figure [TJ] that strong correlations exist between the FP coefficients and the environment 
parameters (-Dec and p), while the correlations are less marked (or absent) with the galaxy 



■^It is the radius at which the mean interior overdensity is 200 times the critical density of the Universe. 

"^The local density around each galaxy has been computed in the circular area containing the 10 nearest 
neighbors with My < —19.5: p = 10/ttRIq {Rio in Mpc). The computation is a bit more complex for the 



objects close to the edge of the WINGS CCD frame s. A statistical background correction of the counts has 



been applied using the recipe bv iBerta et al. 



200^ 
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Fig. 14. — The residuals of the FP fit versus the local density (upper panel) and the 
normalized cluster-centric distance (bottom panel). Note the lack of correlation in this plots 
with respect to that found in Figure [151 
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Fig. 15. — Average values of the the FP coefficients in different bins of color [{B — V)], 
Sersic index [Log(n)], axial ratio [b/a], clustercentric distance [Log{Dcc / R200)] and local 
density [Log(p)] as a function of the median values of the same quantities inside the bins. 
The number of galaxies per bin is fixed to 150 and the BCGs are assumed to coincide with 
the centers of the clusters. The correlation coefficients of the different pairs of bin-averaged 
quantities are also reported in the panels. Note the strong correlations between the FP 
coefficients and the environment parameters Dec and p. 
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properties. We have verified that the average (and median) absolute magnitudes do not 
vary significantly in the different bins of the test quantities Log{Dcc) and Log(p). Thus, 
the correlations among these quantities and the FP coefficients cannot be induced by the 
above mentioned dependence of the FP coefficients on the absolute magnitude (see Sec. I4.2p . 
Moreover, the lack of correlation in the three uppermost panels of Figure [16] rules out the pos- 
sibility that the above trends just reflect similar trends involving the very physical quantities 
that define the FP. 

Figure [I5] suggests that, in the FP, the dependences on both velocity dispersion and 
average surface brightness of galaxies become lower and lower as the distance from the 
cluster center increases. Looking at the two leftmost panels of Figure [15], one could wonder 
if the further dependences of the FP coefficients on both the Sersic index (stronger) and the 
color (weaker) are merely reflecting the correlation with the clustercentric distance. The lack 
of correlation in the two lowest panels of Figure [16] help to clarify this point, suggesting that 
light concentration and color could actually be additional (independent) physical ingredients 
of the FP recipe. It is also interesting to note in Figure [15] that the b coefficient correlates 
quite well with the Sersic index n, while a does not. This is likely because b is the coefficient 
associated with the photometric parameter (/i)e , which is in turn obviously related to the 
concentration index n. Finally, we note that, from the very (linear) expression of the FP, 
most of the dependences of the c coefficient on the various test quantities in Figure [15] are 
likely induced by the corresponding dependences of the a and b coefficients, any increase of 
the last ones producing necessarily a decrease of the former one, and viceversa. 

The trends observed in Figure [T5] further confirm that the FP coefficients depend on the 
particular criteria used in selecting the galaxy sample. They are also likely able to explain the 
large scatter of the FP coefficients which is found even after removal of the luminosity-driven 
bias discussed in Section W?2\ (high values of xl and in Figure [7]). 

5.2. FP versus global cluster properties 

We have also explored the possible dependence of the FP coefficients (in particular 
of the coefficient a) on several measured quantities related to the global cluster properties. 
Tentative correlations have been performed with the velocity dispersion of the galaxies in the 
clusters, with the X-Ray luminosity, with redshift, with the integrated V-band luminosity 
(within My = —19.5) of the clusters, with different kinds of cluster radii, with the average 
Sersic index of the cluster galaxies, with the average Log(M/L), etc. Some of these plots 
are shown in Figure [17] No significant correlations have been found. 
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Fig. 16. — Normalized cluster-centric distance versus velocity dispersions, effective surface 
brightness, effective radii, local density, Sersic index and color for our sample of early-type 
galaxies with My < —19.5. Note the well known (obvious) dependence of the local density on 
the clustercentric distance, while the other panels display the substantial lack of correlations 
with the other variables. The obvious cut at large Log{n) is due to the intrinsic limit of 
GASPHOT to give Sersic index n > 8. 
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Fig. 17. — From top to bottom: the FP coefficient a versus the number of galaxies (with 
My < —19.5) in the cluster, the radius of the cluster /?2oo, the X-ray luminosity, the rms of 
peculiar velocities of galaxies in the cluster, the integrated V-band luminosity (again with 
My < —19.5) and the redshift. No significant correlations are found. 
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5.3. Can we provide a general recipe for deriving the FP ? 

From the analysis performed in Section |H the differences in the FP coefficients appear to 
be related to sampling aspects (i.e. the luminosity cut-off). In Section[5]we have shown that, 
although not depending on the global cluster properties, the FP coefficients are also strongly 
related to the environmental properties of galaxies [Dec and p) and to their internal structure 
(Sersic index). These dependences likely concern the very formation history of galaxies and 
clusters. They are not strictly referable as sampling effects, but we can of course always 
speak of sampling, as far as they translate into the photometric properties of galaxies. This 
let us understand that the various dependences are actually linked each other and it is not 
easy to isolate each of them. Moreover, it is worth stressing that the previous analyses (never 
tried before) have been made possible just because we have at our disposal a huge sample of 
galaxies, obtained putting altogether data from many different clusters. When dealing with 
the determination of the FP for individual (possibly far) clusters, we usually must settle for 
what we actually have, that are a few galaxies (a few dozens, in the most favourable cases) 
with different luminosities and structures, located in a great variety of environments. In this 
cases, we can hardly renounce to each single galaxy and the above dependences turn out to 
be irreparably entangled each one another. 

From the previous remarks it stands to reason that, even with our large sample of galax- 
ies, to provide a general recipe for determining unbiased coefficients of the FP in individual 
clusters is far from being a realistic objective. The best we can do is to remove from our 
global (W+N+S) sample the low-luminosity galaxies (My >-19.5; see Section and to 
provide, for this restricted sample, the FP coefficients obtained with both MIST and ORTH 
fitting algorithms. They are: 

a = 1.097 ± 0.020; b = 0.318 ± 0.004; c = -8.41 ± 0.097 (MIST) 

a = 1.208 ± 0.052; b = 0.318 ± 0.010; c = -8.65 ± 0.19, (ORTH) 

which we assume to define the global (unbiased, as far as possible), V-band FP of early-type 
galaxies in nearby clusters. The MIST coefficients are more suitable for distance determina- 
tion, while the ORTH ones more properly define the physical relation among the quantities 
involved in the FP. 

Columns 10-15 of Table [3] report the FP coefficients of the individual clusters obtained 
running MIST just for galaxies with My <-19.5. 
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THE M/Lv RATIO OF EARLY-TYPE GALAXIES IN NEARBY 

CLUSTERS 



The variation of the mass-to-hght ratio with luminosity is the most popular explanation 
for the tilt of the FP with respect to the virial expectation. Therefore, it is important to 
analyse the M/ L ratio of cluster galaxies with our extensive photometric data, which account 
for the non homologous structure of the ETGs by means of the Sersic parameter n. 



According to ldi Serego et al. I (120051 ) (see also lMichard Ill980l ). we calculate the dynam- 
ical mass of galaxies using the formula: M/Mq = Kv(n) a^Rp/G, where the v irial coefficient 



2OO2I) and G is the 



Kyln) is a decreasing function of the Sersic index n (iBertin et al. 
gravitational constant. 

Figure [18] shows the mass-to-light ratio as a function of mass for our global galaxy 
sample. The full straight line in the figure represents the linear fit obtained minimizing the 
weighted sum of perperdicular distances from the line itself [see the equa tion F(M) in the 



figure] . The open dots refer to the sample of galaxies in Coma given by |j0rgensen et al. 



(Il996l ). with the relative orthogonal fit represented by the dashed line. Note in particular 
that in our global sample the scatter of the residuals relative to the best-fit relation is greater 
than in the Coma sample (0.19 vs. 0.11; for the NFF S and SPSS samp l es the scatters are 



0.18 and 0.21, respectively). However, we recall that iJeirgensen et al. I ( 1l996l ) derived the 
photometric parameters {Re and {fi)e ) and the mass by assuming r^/^ luminosity profiles, 
while we adopted the more general Sersic profiles. This might also explain the fact that 
the slope of the relation for our global sample [0.511(±0.019)] is much larger that in the 
Coma sample [0.28(±0.028)]. By the way, the slopes we found for the NFPS and SDSS 
samples separately, are quite consistent each other, within the errors [0.522(±0.022) and 
0.600(±0.052), respectively]. 



Figure [19] reports several plots showing the correlations among different measured and 
evaluated quantities involving the mass-to-light ratio estimate. In particular we test (at 
the ordinate) dynamical masses, mass-to-light ratios and residuals [Log(M/L)-F(M)] of the 
relation in Figure [18] versus (at the abscissa) Sersic indices, velocity dispersions, effective 
radii and luminosities. 

Some of the correlations in Figure [TH] are well known (i.e. mass vs. luminosity), obvious 
(i.e. M/L residuals vs. luminosity), or expected by definition (i.e. mass vs. a and -Rc )• 
Less obvious seem to be some other correlations (i.e. M/L vs. cr and M/L residuals vs. 
a and Re ) or lack of correlations (i.e. M vs. Sersic index, M/L vs. Re and luminosity). 
For instance, according to the formula we used to derive the dynamica l mass, it should be a 
strongly decreasing function of the Sersic index (see lBertin et al. II2002I ). while the correlation 
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Fig. 18. — Mass-to-light ratio versus dynamical mass for our global galaxy sample. Black 
and gray (green in the electronic version) dots refer to W-l-N and W+S galaxies, respectively, 
while open dots refer to a sample of galaxies in Coma (see text). The solid line gives the 
orthogonal fit of the W-l-N+S data, while the dashed line reports the fit for the Coma sample. 
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Fig. 19. — Correlations among measured (abscissa) and evaluated (ordinate) quantities 
involved in the mass-to-light ratio estimate. In each panel the proper correlation coefficient 
is also reported. Symbols are as in Figure [181 The straight line reported in the plot Log(M) 
- Logi?e represents the orthogonal best-ffi of the data we discuss in the last paragraph of 
this Section. 
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coefficient of the plot M-n in Figure [19] is slightly positive. Moreover, in the same figure 
the M/L ratio does not seem to correlate at all with either radius or luminosity, while a 
correlation M / L-L has been often claimed to explain the 'tilt' of the FP. 

In Figure [19] we find of particular interest the correlation between M/L and Sersic index 
and that between M/L residuals and Sersic index. The first one, coupled with the lack of 
correlation between M/ L and luminos i ty (w hich is indeed expected for non homologous 



ETGs, as suggested by iTrujillo et al. 1 12004| ). indicates that, for a given luminosity, the 



galaxies showing lower light concentration (lower Sersic index) are more massive (more dark 
matter?). 

The second, even stronger correlation is quite interesting as well. In fact, from the very 
definition of the M/L residuals of the relation in Figure [TSl for a given dynamical mass, the 
lower the residual, the brighter the galaxy. Therefore, the correlation in Figure [TH] between 
M/L residuals and Sersic indices implies that (again for a given mass) the higher the light 
concentration (Sersic index), the brighter the galaxy. 

Thus, the picture emerging about the infiuence of the light concentration in determining 
dynamical mass and luminosity of ETGs is that: (i) for a given luminosity, the higher the 
light concentration, the lower the dynamical mass; (ii) for a given dynamical mass, the higher 
the light concentration, the higher the luminosity. This twofold dependence on the Sersic 
concentration index is expressed by the linear equation: 

Log(n)=1.60xLog(L)-1.16xLog(M)-2.93, 

we have derived minimizing the orthogonal distances from the fitting plane of the points 
in the parameter space {n,L,M). Note that the correlation coefficient between the Sersic 
index computed from this equation and the measured one is CC= 0.59. 

Still concerning the infiuence of the light concentration on the mass-to-light ratio of 
early-ty pe galaxies, it is well known that the Sersic index n correlates with the velocity dis- 



persion (iGraham Il2002[ ). Thus, it is not meaningless wondering if the correlations involving 
n in Figure [T^] just refiect the correlations with a. The upper and middle panels of the 
figure clearly rule out this hypothesis as far as the correlations with mass and mass-to-light 
ratio are concerned (both are positive for a, while for n they are close to zero and negative, 
respectively). Instead, the correlations of the M/L residuals with n and a (lower panels) 
have the same sign. It is worth noting, however, that the correlation turns out to be tighter 
with n than with a and that the same happens (even if with opposite trends) also for the 
M/L ratio in the middle panels of the figure. This might suggest that the driving parameter 
for M/L is actually the light concentration and that the trends with a are just consequence 
of that. 
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We have previously guessed that the different slopes we find in the relation [M/L-M) 
between our sample and the Coma sample could be at least partially due to the different 
models of luminosity profiles used to derive the photometric parameters of galaxies (Sersic 
law and r^/^ law, respectively). Now, we could legitimately guess that the correlations 
shown in the leftmost panels of Figure [19] are artificially produced by the use of the Sersic 
law in deriving the parameters Re and Ky, involved in the computation of the galaxy mass- 



Actu ally, Ky turns out to be a decreasing function of the Sersic index (see iBertin et al. 



2OO2I ). just like the M/L ratio and the M/L residuals in Figure [H] (but, in the same figure 



note the direct, although weak, correlation between the mass and the Sersic index!). 

Trying to clarify these points, we have re-calculated the masses and the luminosities 
of the early-type galaxies in the original W+S sample (397 objects, before selection on a 
and My) using the surface photometry parameters provided by the SDSS database {r^^^ 
profiles). In Figure [20] we plot the M/L ratio and the residuals of the M/L-M relation 
versus the Sersic index for the W+S galaxy sample alone. Left and right panels illustrate 
the relations obtained when masses and luminosities are computed using the Sersic and r^/^ 
surface photometry parameters, respectively. In the former case, the correlation coefficients 
turn out to be undistinguishable from those obtained for the whole (W+N+S) galaxy sample 
(see Figure [T9]) . In the latter case the correlations are less pronounced, but still they are in 
place, as indicated by 10,000 random extractions of couples of uncorrelated vectors having 
the same dimension (397) and distributions of the real ones. In fact, the probability that 
the correlation coefficients of the real sample in the right panels of Figure [20] are drawn from 
a parent population of uncorrelated quantities turns out to be very small: ~0.005 and ~0 
for the correlations in the upper-right and lower-right panels, respectively. This enforces our 
previous conclusions about the dependence of masses and luminosities of early-type galaxies 
on the Sersic index. The weaker correlations found with the r^/^ profiles if compared with 
the Sersic profiles, are likely the consequence of having forced the real luminosity structure 
of galaxies to obey the de Vaucouleurs law. Finally, we mention that the slope of the relation 
(M/L-M) for the W+S sample turns out to be 0.47 and 0.38 with the Sersic- and r^/^- 
law approaches, respectively, thus confirming our previous guess that it is influenced by the 
assumption about the luminosity profile of galaxies (see the comparison between our sample 
and the Coma sample in Figure [T5I) . 



In a recent paper [Robertson et al. I (120061 ) claim that the tilt of the FP is closely linked 



to dissipation effects during galaxy formation. They show the results of their simulations 
in a plot of M/Mq vs. the radius {R) of the galaxies. Figure \W\ also shows a similar plot 
for our galaxy sample. Although a correlation between M and Re is expected from the very 
definition of dynamical mass, we note in this plot that the more massive objects are preferably 
found below the orthogonal best-fit of the data distribution (by the way, these objects are 
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Fig. 20. — The M/ L ratio and the residuals of the M/ L-M relation versus the Sersic index 
for the W+S galaxy sample alone. Left and right panels illustrates the relations obtained 
using the Sersic and i?^/^ surface photometry parameters, respectively, in the computation 
of masses and luminosities. Again, the correlation coefficients are reported in each panel. 
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those deviating in the high luminosity region from the FP projection in Figure [T3l) . High 
luminosity objects also deviate with respect to the bulk of the early-type p opulation in the 



R — L relation, that for our global sample has the same slope found by iBernardi et al. 



( 120071 ) (i? oc L"-^^). The systematically larger size of these galaxies in th at relation may 
be consistent with the results of the simulations by [Robertson et al. I (120061 ). if one invokes 
the dry dissipationless merger mechanism for their formation. It again points towards the 
hypothesis that the FP relation might be non linear, this time in its high-mass region (see 
Section U2] for a similar finding in the low-mass region). Thus, in the parameter space of the 
FP the real surface defined by ETGs could be slightly bent, reflecting the different formation 
mechanisms producing the present day ETGs. In Section 14.21 we have seen that part of the 
scatter of the FP coefficients is just due to such an effect, coupled with the different statistical 
properties of the galaxy samples. As a consequence, to draw any conclusion about luminosity 
evolution and downsizing effect might be dangerous (the slopes of the relations depend on 
the sample selection rules/biases), unless the galaxy samples involved in these analyses, even 
spanning wide ranges of redshift, share the same distributions of photometric/kinematic 
properties of galaxies. 



7. CONCLUSIONS 

We have derived the Fundamental Plane of early-type galaxies in 59 nearby clusters 
(0.04 < z < 0.07) by exploiting the data coming from three big surveys: WINGS(W), to 
derive Re and (/i)e , and NFPS+SDSS(N-|-S), to derive a. The fits of the FP, obtained 
for the global samples W+N and W+S, as well as for each cluster, have revealed that the 
FP coefficients span considerable intervals. By means of extensive simulations, we have 
demonstrated that this spread is just marginally consistent with the statistical noise due to 
the limited number of galaxies in each cluster. It seems at least partly due to a luminosity- 
driven bias depending on the statistical properties of the galaxy samples. These can be 
induced both by observing limitations and by selection rules. In fact, even if the best-fitting 
solution obtained for the global W+N+S dataset does not differ significantly from previous 
determinations in the literature, systematic different FP coefficients are found when the 
galaxy samples are truncated in the faint-end part at different cut-off absolute magnitudes. 
We speculate that, rather than a plane, the so called FP is actually a curved surface, which 
is approximated by different planes depending on the different regions of the FP space 
occupied by the galaxy samples under analysis. To this concern, we could go farther on in 
the speculation, suggesting that a bent FP coul d be, at least partially , reconciled with the 



numerical simulations in ACDM cosmology (see iBorriello et al. II2003I ). By the way, such a 



speculation could also be supported by the large scatter we find in the M/ L-M relation, at 
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variance with other determinations, whose tightness has been sometime invoked to rule out 
the hierarchical scenario. 

Perhaps the most interesting result of the present analysis concerns the dependence of 
the FP coefficients on the local environment, which clearly emerges when we derive the FP 
in different bins of the cluster-centric distance and local density. Finally, we do not find any 
dependence of the FP coefficients on the global properties of clusters. 

Concerning the M/L ratio, we also find that both M/L and the residuals of the M/L- 
M relation turn out to be anti-correlated with the Sersic indices. These trends could imply 
that, for a given luminosity, more massive galaxies display a lower light concentration, while 
for a given dynamical mass, the higher the light concentration, the brighter the galaxy. 

The main results of this work can be summarized as follows: 

• the FP coefficients depend on the adopted fitting technique and (marginally) on the 
methods used to derive the photometric parameters and (/x)e ; 

• The observed scatter in the FP coefficients cannot be entirely ascribed to the uncer- 
tainties due to the small number statistics; 

• the FP coefficients depend on the distributions of photometric/kinematic properties of 
the galaxies in the samples (mainly on the faint-end luminosity cut-off); 

• the FP coefficients are strongly correlated with the environment (cluster-centric dis- 
tance and local density), while the correlations are less marked (or absent) with the 
galaxy properties (Sersic index, color and flattening); 

• the FP coefficients do not correlate with the global properties of clusters (radius, ve- 
locity dispersion. X-ray emission, etc.); 

• the distribution of galaxies in the FP parameter space suggest that the variables R^, 
(/x)e , and a define a slightly warped surface. Forcing this surface to be locally a plane 
causes a systematic variation the FP coefficients, depending on the selection rules used 
to define the galaxy sample; 

• using the FP as a tool to derive the luminosity /size evolution of ETGs may be dan- 
gerous, unless the galaxy samples involved in the analysis are highly homogeneous in 
their average photometric properties; 

• the M/L ratio is not correlated with L when the non homology of ETGs is taken 
into account. This is an indication that most of the tilt of the FP is indeed due to 
dynamical and structural non-homology of ETGs; 
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• the mutual correlations among mass, luminosity and light concentration of ETGs indi- 
cate that, for a given mass, the greater the light concentration the higher the luminosity, 
while, for a given luminosity, the lower the light concentration, the greater the mass; 

• the bending of the FP and the large scatter found in the M/L-M relation could, at 
least partially, reconcile the FP phenomenology with the hierachical merging scenario 
of galaxy formation. 



By exploiting the galaxy mass estimates coming from both the K-ban d WINGS data an d 



the spectro-photometric analysis of the galaxies in the WINGS survey (IFritz et al. 1120071 ) 



in a following paper we will go into more depth about the scaling relations involving mass, 
structure and morphology of galaxies in nearby clusters. 



We wish to thank the anonymous referee, since her/his useful comments helped us to 
improve the final version of the paper. We also wish to thank our colleagues from the 
Astronomical Observatory of Napoli (La Barbera et al. 2000) for having kindly provided us 
with the fitting tool MIST, which we have intensively used in this paper to derive the FP 
coefficients. 
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Table 1. The Cluster sample. 
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A1795 


0.0625 


883 


45.05 


-25.57 


1.978 


-23.56 


27 








A1831 


0.0615 


565 


44.28 


-25.63 


1.266 


-22.93 


21 


_ 


_ 


A1983 


0.0436 


563 


43.67 


-24.87 


1.272 


-22.08 


14 


_ 


_ 


A1991 


0.0587 


557 


44.13 


-25.51 


1.250 


-23.23 


20 


_ 


_ 


A2107 


0.0411 


634 


44.04 


-24.90 


1.435 


-23.28 


27 






A2124 


0.0656 


885 


44.13 


-25.51 


1.980 


-23.53 


30 


38 


19 


A2149 


0.0650 


393 


43.92 


-25.61 


0.879 


-23.24 




20 




A2169 


0.0586 


529 


43.65 


-24.80 


1.188 


-22.49 




10 




A2256 


0.0581 


1353 


44.85 


-26.37 


3.038 


-23.40 


33 






A2382 


0.0618 


998 


43.96 


-25.75 


2.234 


-22.84 


20 






A2399 


0.0579 


781 


44 


-25.32 


1.754 


-22.60 


24 


25 


15 


A2572a 


0.0403 


650 


44.01 


-24.94 


1.472 


-23.26 


10 






A2589 


0.0414 


972 


44.27 


-24.78 


2.200 


-23.45 


22 






A2593 


0.0413 


729 


44.06 


-24.97 


1.650 


-22.84 




23 




A2657 


0.0402 


673 


44.2 


-24.85 


1.524 


-22.68 


21 
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Table 1 — Continued 



Cluster 


z(NED) 


S 


Lx 


Mtot 


^200 


Mv{BCG) 


Nw+N 


Nw+s 


^W+[NgcS] 


A 07'iA 


U.UD/O 


sn/i 


44. 4i 


-ZO.UD 


1 sno 

l.oUZ 


OQ /IS 
-Zo. 40 


OS 
Zo 






A Q1 OS 
A.0LZ0 


u.uoyy 


y ( D 


AA 


Of? Ofi 
-ZD. ZD 


1 on 
z.iyu 


OQ 
-Zo.OO 


A7 
4 ( 






A 01 CO 
AoiOo 


u.uoy ( 


ill ( 


/I /I 7"? 
44. ( 


OR 1 Q 
-ZD.iO 


i^n7 

Z.OU ( 


OQ so 

-Zo.oZ 


A 1 
41 






AoZDO 


U.Uooy 


1400 


A A 70 

44. ( y 


Of? OQ 

-ZD.Zo 


Q OQQ 


OQ CO 

-Zo.oy 


A n 
4U 






A QQ7fi 


n n/i 

U.U40D 


yuz 


AA 

44. oy 


01^ fl/l 
-Z0.U4 


n'^7 

Z.Uo ( 


0'^ 1 

-Zo. iz 


on 
zu 






A ^"^Qc; 


U.UOUD 


1 1 Qc; 


AA Af^ 
44.40 


oc; 07 
-zo.y 1 


ROO 

z.oyz 


-zo.oy 


o4 






A ^/IQ7 


U.UD / 1 


787 


AA ^fi 
44. ID 


oc; c;^ 

-ZO.Oo 


1 7c;q 
1. 1 oy 


00 /I c; 

-ZZ. 40 


1 R 
±D 








U.UOoO 


iuyo 


/I /I 1 9 
44. IZ 


0^ 71 
-ZO. ( i 


/I P^Q 

z.4oy 


0"? 77 
-Zo. / ( 


OQ 
Zo 






AoOZOD 


U.Uooo 


070 

y ( y 


44.0 


op; 7/1 
-ZO. ( 4 


onQ 
Z.ZUo 


OQ fil 

-Zo. Di 


1 /I 
14 






A Qc^Qn 
AoOoU 


U.Uuo 1 


DOO 


A Q 0/1 

4o.y4 


op; p;p; 
-ZO.Oo 


1 p^/l 1 


OQ 7Q 

-Zo. / 


Of? 
ZD 






A QC^QO 


U.U004 


7p;n 
1 oU 


A A AK 
44.40 


OK 01 

-zo.yi 


i.DOD 


OQ 7n 
-Zo. / U 


Q7 
f 






A QKCCf! 
AoOOD 


u.U4< y 


044 


A'i 07 
4o.y 1 


oc; c;s 

-ZO.Oo 


1 /I 

1.40O 


OQ Q/l 
-Zo.o4 


oc; 
zo 






A3558 


0.0480 


989 


44.8 


-26.38 


2.232 


-24.18 


52 






A3560 


0.0489 


844 


44.12 


-25.68 


1.903 


-22.09 


19 






A3667 


0.0556 


1170 


44.94 


-26.15 


2.631 


-23.97 


54 






A3716 


0.0462 


855 


44 


-25.95 


1.932 


-22.94 


37 






A3809 


0.0620 


631 


44.35 


-25.35 


1.414 


-22.85 


27 






A3880 


0.0584 


893 


44.27 


-25.02 


2.005 


-23.07 


16 






A4059 


0.0475 


843 


44.49 


-25.25 


1.901 


-23.64 


29 






IIZW108 


0.0493 


579 


44.34 


-25.41 


1.306 


-23.77 


30 






MKW3s 


0.0450 


575 


44.43 


-24.69 


1.299 


-22.72 


22 






RX1022 


0.0534 


777 


43.54 


-25.16 


1.748 


-22.66 




11 




RX1740 


0.0430 


596 


43.7 


-24.27 


1.347 


-22.41 


11 






Z2844 


0.0500 


559 


43.76 


-23.93 


1.260 


-23.31 


21 






Z8338 


0.0473 


747 


43.9 


-25.06 


1.684 


-23.15 


12 






Z8852 


0.0400 


795 


43.97 


-25.30 


1.800 


-23.41 


18 
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Table 2: FP coeff. for different galaxy samples, fitting algorithms and surface photometries. 



Sample 


a 


b 


c 


rmsa 


rmsb 


rmsc 




Fitting 


Phot. 


W+N+S 
W+N+S 


1.152 
1.293 


0.320 
0.322 


-8.56 
-8.91 


0.021 
0.021 


0.004 
0.003 


0.095 
0.002 


1550 
1550 


MIST 
ORTH 


WINGS 
WINGS 


W+N 
W+N 


1.113 
1.258 


0.319 
0.329 


-8.45 
-8.99 


0.021 
0.022 


0.004 
0.003 


0.102 
0.003 


1368 
1368 


MIST 
ORTH 


WINGS 
WINGS 


W+S 
W+S 


1.332 
1.306 


0.318 
0.303 


-8.93 
-8.56 


0.050 
0.048 


0.008 
0.008 


0.198 
0.006 


282 
282 


MIST 
ORTH 


WINGS 
WINGS 


W+S 


1.297 


0.319 


-8.87 


0.050 


0.008 


0.198 


282 


MIST 


SDSS 


COMA 
COMA 


1.239 
1.439 


0.342 
0.345 


-9.15 
-9.67 


0.080 
0.077 


0.013 
0.013 


0.310 
0.010 


80 
80 


MIST 
ORTH 


JORG 
JORG 



Table 3. FP MIST coefficients of the individual clusters. 







All Galaxies 








Galaxies with My <-19.5 




Cluster Ng 


a 


rms(a) b rms(b) 


c 


rms(c) Ng 


a 


rms(a) b rms(b) c 


rms(c) 



A0085 


63 


1 


137 





083 





304 


A1069 


20 


1 


236 





216 





275 


A119 


48 


1 


289 





127 





289 


A1291 


13 


1 


415 





222 





381 


A133 


23 


1 


462 





161 





371 


A160 


23 


1 


580 





247 





350 


A1631a 


22 


1 


214 





093 





289 


A 1644 


41 


1 


030 





114 





323 


A1668 


23 





781 





192 





274 


A168 


12 


1 


279 





129 





344 


A 1795 


27 





774 





120 





255 


A1831 


21 





725 





116 





325 


A 1 QCQ 


1 A 


1 
i 




u 


1 f^R 
100 


u 


071 


A1991 


20 





933 





183 





253 


A2107 


27 





993 





127 





269 


A2124 


49 


1 


065 





131 





317 


A2149 


20 


1 


159 





127 





321 


A2169 


10 


1 


500 





165 





331 


A2256 


33 





885 





129 





302 


A2382 


20 


1 


547 





199 





314 


A2399 


34 


1 


154 





113 





338 


A2572a 


10 


1 


157 





267 





299 


A2589 


22 


1 


016 





188 





315 


A2593 


23 


1 


559 





221 





320 



0.013 


-8.24 


0.33 


52 


1.013 


0.030 


-7.81 


0.75 


20 


1.236 


0.018 


-8.21 


0.54 


45 


1.169 


0.016 


-10.37 


0.62 


11 


1.635 


0.009 


-10.27 


0.41 


23 


1.462 


0.040 


-10.10 


0.89 


19 


1.743 


0.014 


-8.06 


0.23 


21 


1.187 


0.019 


-8.35 


0.56 


40 


1.088 


0.030 


-6.80 


0.90 


23 


0.781 


0.058 


-9.27 


1.16 


11 


1.111 


0.029 


-6.41 


0.79 


27 


0.774 


0.020 


-7.74 


0.38 


21 


0.725 


0.039 


-7.33 


0.83 


13 


0.995 


0.051 


-6.68 


1.19 


20 


0.933 


0.028 


-7.22 


0.64 


25 


0.981 


0.014 


-8.27 


0.50 


48 


0.992 


0.023 


-8.58 


0.62 


20 


1.159 


0.047 


-9.55 


0.78 


8 


1.441 


0.018 


-7.55 


0.43 


33 


0.885 


0.019 


-9.32 


0.54 


20 


1.547 


0.017 


-8.91 


0.39 


30 


1.002 


0.047 


-8.20 


1.45 


10 


1.157 


0.038 


-8.20 


1.06 


20 


0.854 


0.052 


-9.48 


1.28 


15 


0.698 



0.076 


0.289 


0.014 


-7.64 


0.36 


0.216 


0.275 


0.030 


-7.81 


0.75 


0.113 


0.287 


0.018 


-7.90 


0.52 


0.389 


0.377 


0.018 


-10.78 


1.09 


0.161 


0.371 


0.009 


-10.27 


0.41 


0.303 


0.363 


0.050 


-10.73 


1.22 


0.106 


0.290 


0.014 


-8.02 


0.23 


0.124 


32Q 


020 


-8 61 


60 


0.192 


0.274 


0.030 


-6.80 


0.90 


0.092 


0.316 


0.051 


-8.32 


1.03 


0.120 


0.255 


0.029 


-6.41 


0.79 




n 325 


020 


-7 74 


38 


U.ioU 


yj.iot 


n c\Ar\ 
U.U4U 


-D.yo 


U.oo 


0.183 


0.253 


0.051 


-6.68 


1.19 


0.151 


0.296 


0.014 


-7.72 


0.46 


0.114 


0.314 


0.014 


-8.05 


0.44 


0.127 


0.321 


0.023 


-8.58 


0.63 


0.136 


0.286 


0.053 


-8.48 


0.98 


0.129 


0.302 


0.018 


-7.55 


0.44 


0.199 


0.314 


0.019 


-9.32 


0.54 


0.092 


0.324 


0.014 


-8.28 


0.33 


0.267 


0.299 


0.047 


-8.20 


1.45 


0.163 


0.303 


0.039 


-7.59 


1.07 


0.089 


0.189 


0.028 


-4.91 


0.74 



Table 3 — Continued 







All Galaxies 








Galaxies with My <-19.5 




Cluster Ng 


a 


rms(a) b rms(b) 


c 


rms(c) Ng 


a 


rms(a) b rms(b) c 


rms(c) 



A2657 


21 


1 


059 





165 





336 


A2734 


28 


1 


071 





198 





325 


A3128 


47 


1 


329 





133 





365 


A3158 


41 


1 


189 





090 





306 


A3266 


40 





976 





105 





337 


A3376 


20 


1 


174 





210 





293 


A3395 


34 


1 


066 





091 





374 


A3497 


16 





731 





151 





274 


A3528a 


23 





730 





150 





334 


A3528b 


14 


1 


246 





163 





256 


A3530 


26 





956 





107 





313 


A3532 


37 


1 


101 





099 





326 




oc; 
zo 


1 
i 


00/1 


u 


1 71 


u 




A3558 


52 


1 


014 





079 





360 


A3560 


19 


1 


284 





225 





303 


A3667 


54 


1 


208 





089 





326 


A3716 


37 


1 


277 





166 





323 


A376 


27 


1 


067 





108 





307 


A3809 


27 





903 





073 





329 


A3880 


16 


1 


096 





115 





397 


A4059 


29 


1 


149 





134 





336 


A548b 


22 





991 





120 





325 


A602 


18 


1 


180 





231 





361 


A671 


16 


1 


123 





161 





273 



0.031 


-8.70 


0.84 


21 


1.059 


0.028 


-8.52 


0.93 


28 


1.071 


0.024 


-9.86 


0.63 


47 


1.329 


0.019 


-8.34 


0.40 


41 


1.189 


0.015 


-8.51 


0.40 


40 


0.976 


0.032 


-8.07 


1.01 


20 


1.174 


0.021 


-9.47 


0.44 


34 


1.066 


0.018 


-6.64 


0.59 


16 


0.731 


0.019 


-7.90 


0.61 


23 


0.730 


0.021 


-7.53 


0.64 


13 


1.236 


0.017 


-8.03 


0.35 


26 


0.956 


0.020 


-8.60 


0.42 


37 


1.101 


0.028 


-10.08 


0.76 


24 


1.155 


0.016 


-9.06 


0.38 


52 


1.014 


0.046 


-8.52 


0.78 


19 


1.284 


0.019 


-8.82 


0.47 


54 


1.208 


0.019 


-8.85 


0.48 


37 


1.277 


0.019 


-8.03 


0.51 


27 


1.067 


0.017 


-8.18 


0.33 


27 


0.903 


0.057 


-9.96 


1.22 


16 


1.096 


0.020 


-8.91 


0.58 


26 


1.127 


0.015 


-8.36 


0.52 


18 


0.941 


0.048 


-9.40 


1.37 


16 


0.948 


0.020 


-7.58 


0.64 


14 


1.175 



165 


336 


031 


-8 70 


n 84 


0.198 


0.325 


0.028 


-8.52 


0.93 


0.133 


0.365 


0.024 


-9.86 


0.63 


0.090 


0.306 


0.019 


-8.34 


0.40 


0.105 


0.337 


0.015 


-8.51 


0.40 


0.210 


0.293 


0.032 


-8.07 


1.01 


0.091 


0.374 


0.021 


-9.47 


0.44 


0.151 


0.274 


0.018 


-6.64 


0.59 


0.150 


0.334 


0.019 


-7.90 


0.61 


0.197 


0.256 


0.025 


-7.49 


0.84 


0.107 


0.313 


0.017 


-8.03 


0.35 


099 


326 


020 


-8 60 


0.42 


U.i ( t 


U.ooi 


n noT 


-y.oo 


U. (o 


0.079 


0.360 


0.016 


-9.06 


0.38 


0.225 


0.303 


0.046 


-8.52 


0.78 


0.089 


0.326 


0.019 


-8.82 


0.48 


0.166 


0.323 


0.019 


-8.85 


0.48 


0.108 


0.307 


0.019 


-8.03 


0.51 


0.073 


0.329 


0.017 


-8.18 


0.33 


0.115 


0.397 


0.057 


-9.96 


1.22 


0.137 


0.343 


0.024 


-9.00 


0.66 


0.082 


0.317 


0.010 


-8.09 


0.32 


0.175 


0.330 


0.038 


-8.25 


1.02 


0.218 


0.296 


0.028 


-8.16 


0.96 



Table 3 — Continued 











All Galaxies 










Galaxies with My <-19.5 






Cluster 




a 


rms(a) 


b 


rms(b) 


c 


rms(c) 




a 


rms(a) 


b 


rms(b) 


c 


rms(c) 


A754 


46 


0.993 


0.107 




U.UID 


-8.17 


0.40 


46 


0.993 


n 1 fi7 
u.iu / 


n ^1 7 


U.UID 


-8.17 


0.40 


A780 


16 


1.364 


0.215 


0.325 


0.032 


-9.14 


0.96 


16 


1.364 


0.215 


0.325 


0.032 


-9.14 


0.96 


A957x 


29 


1.266 


0.113 


0.319 


0.014 


-8.80 


0.39 


23 


1.085 


0.093 


0.312 


0.012 


-8.24 


0.31 


A970 


25 


1.156 


0.223 


0.317 


0.022 


-8.49 


0.73 


25 


1.156 


0.223 


0.317 


0.022 


-8.49 


0.72 


IIZW108 


30 


0.957 


0.137 


0.244 


0.031 


-6.64 


0.79 


29 


0.932 


0.138 


0.243 


0.031 


-6.57 


0.78 


MKW3s 


22 


1.112 


0.151 


0.260 


0.024 


-7.31 


0.59 


20 


1.099 


0.132 


0.259 


0.028 


-7.24 


0.62 


RX1022 


11 


1.230 


0.146 


0.305 


0.030 


-8.47 


0.30 


9 


1.138 


0.147 


0.291 


0.056 


-7.97 


1.04 


RX1740 


11 


1.327 


0.408 


0.427 


0.047 


-11.08 


1.41 


11 


1.327 


0.408 


0.427 


0.047 


-11.08 


1.41 


Z2844 


21 


0.961 


0.192 


0.213 


0.021 


-6.00 


0.54 


18 


0.935 


0.201 


0.226 


0.022 


-6.21 


0.53 


Z8338 


12 


0.933 


0.244 


0.273 


0.047 


-7.12 


1.27 


10 


0.685 


0.194 


0.230 


0.033 


-5.69 


0.98 


Z8852 


18 


0.761 


0.097 


0.339 


0.023 


-8.13 


0.52 


17 


0.731 


0.115 


0.337 


0.023 


-8.02 


0.54 
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Table 4: Statistics of the measured coefficients for the MIST ffis of the FP. 



Samnle 


doefficient 


A vera f?e 


St. Dev. 


IVTediaTi 


Notes 




<x 


1 1 21+0 027 


n 297 


1 123 


(all plnci'l'prc; 1 




b 


0.316±0.005 


0.040 


0.319 






c 


-8.41±0.137 


1.052 


-8.35 




W+N+S 


a 


1.108±0.037 


0.207 


1.071 


(clusters with Ng > Nmed) 




b 


0.321±0.006 


0.033 


0.323 






c 


-8.49±0.172 


0.956 


-8.49 




W+N 


a 


1.081±0.029 


0.211 


1.066 


(all clusters) 




b 


0.311±0.006 


0.041 


0.315 






c 


-8.23±0.140 


1.012 


-8.20 




W+N 


a 


1.047±0.033 


0.176 


1.064 


(clusters with Ng > N^^d) 




b 


0.319±0.006 


0.034 


0.323 






c 


-8.31±0.180 


0.953 


-8.34 




W+S 


a 


1.308±0.052 


0.195 


1.279 


(all clusters) 




b 


0.327±0.008 


0.030 


0.324 






c 


-9.01±0.252 


0.913 


-8.73 




W+S 


a 


1.226±0.115 


0.230 


1.201 


(clusters with Ng > Nmed) 




b 


0.313±0.005 


0.011 


0.320 






c 


-8.59±0.335 


0.671 


-8.08 





